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ABSTRACT 

Context. The ephemeris of Phoebe, the ninth satellite of Saturn, is not very accurate. Previous dynamical models were usually too 
simplified, the astrometry is heterogeneous and, the Saturn's ephemeris itself is an additionnal source of error. 

Aims. The aim is to improve Phoebe's ephemeris by using a large set of observations, correcting some systematic errors and updating 
the dynamical model. 

Methods. The dynamical model makes use of the most recent ephemeris of planets and Saturnian satellites. The astrometry of Phoebe 
is improved by using a compilation of ground-based and space-based observations and by correcting the bias in stellar catalogues 
used for the reduction. 

Results. We present an accurate ephemeris of Phoebe with residuals of 0.45 arcsec and with an estimated accuracy of Phoebe's 
position of less that 100 km on 1990-2020 period. 
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1. Introduction 

Phoebe is the biggest irregular satellite of Saturn (220km in di- 
ameter). It has a retrograde orbit with an inclination of about 
176 degrees on the eclipitic plane orbiting at about 13 million 
km (0.086 au) from Saturn. The ephemeris of Phoebe is not very 
accurate for three reasons. The main reason is the astrometry of 
Phoebe. Because of its faintness (about magnitude 16), Phoebe 
is not easy to observe and was especially difficult to observe in 
the past. Most of the available observations have been realized 
in the last twenty years. There are only a few observations in 
the past and they are not very accurate. Another reason is that 
the dynamical models of Phoebe's motion are usually not com- 
plete and take into account only a few perturbations. However, 
while the astrometry is not good, a simple dynamical model 
could be enough to compute the position of Phoebe. The last 
reason is the dependance on Saturn's position. When comparing 
the observed and computed positions of Phoebe, the position of 
Saturn, or more exactly the barycentre of the Saturnian system 
is required because the dynamical model provides the position 
of Phoebe in relation to the Saturnian barycentre. Consequently, 
the ephemeris of Saturn can be the source of systematic error in 
Phoebe's position. 

In this paper, we develop a new ephemeris of Phoebe by im- 
proving the three previous limiting factors. In particular, we first 
develop a new flexible dynamical model that can take into ac- 
count the perturbations of the planets, the main Saturnian satel- 
lites, and the flatness of Saturn, using the most recent planetary 
and satellite ephemerides (Sect.|2|. 



An important step is realized for astrometry. We compile 
a large part of ground-based observations and include also 
observations from Voyager and Cassini spacecrafts (Sect. [3j. 
Moreover, we study old observations and try to reduce some 
old observations with modern techniques (Sect. 4. 1 1. In addi- 



tion, most of the observations (78%) are also corrected of bias in 
the stellar catalogue used for their reduction (Sect. [42]). 

After presenting the fitting process (Sect. [5]), we compare 
four different models of Phoebe's motion (Sect. [6} and high- 
light the one with the best residuals (Sect. [71. Finally, by using 
statistical methods, we also compute the estimated accuracy of 
Phoebe's position during the 1875-2020 period (Sect. [8}. 



2. A flexible dynamical model 



According to |Jacobson] ( |1998| l, the previous dynamical mod- 
els usually took into account the perturbation of some plan- 
1979) includes only the Sun's perturbation, Bykova 



ets. [Rose 
|& Shikha 



ev| (|1982|) ad ded Jupiter's perturbation and |Bec 
Borsenberger & Rocher ( 1982[ ) also added the perturbation of 
Titan. For more recent models, |Jacobson] ( |1998| ) took into ac- 
count the perturbations of the Sun, Jupiter, and Uranus, as well as 
the perturbations of the main Saturnian satellites for post-1966 
observations and the perturbation of an orbiting Titan for pre- 
1966 observations. Shen et al. (2005 2011) included the pertur- 
bations of the Sun, Jupiter and Uranus and the flatness (J2) of 



Saturn, whereas Emelyanov (2007) took into account the pertur- 
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bations of the Sun, Jupiter, Uranus, Neptune and Saturn's flat- 
ness (J2, Ja). For the two last models, the Saturnian satellites are 
taken into account by adding their masses to the mass of Saturn 
and by correcting the J2 and J4 values for Em elyanov] ( |2007| ). 
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Finally, |Jacobson et al.| ( |2006| l proposed a complete model of 
the Saturnian system by determining the motion of all Saturnian 
satellites including Phoebe. 

Figure [T] represents the ratio of maximum and minimum 
magnitude of gravitational accelerations compared to Saturn's 
acceleration. The main perturbers of Phoebe's motion are the 
Sun, Jupiter, Titan, Uranus, and Iapetus. Taking into account the 
flatness of Saturn does not seem justified since other perturba- 
tions like those due to the main Saturnian satellites are not in- 
cluded. In particular, the perturbation due to J4 is less important 
than the gravitational perturbations of the main asteroids. 



1 

0.01 
0.0001 

1e-06 
1e-08 
1e-10 
1e-12 
1e-14 
1e-16 
1e-18 



Table 1. Dynamical constants of Saturn's system ( Jacobson et al. 

|2006| ) 
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Fig. 1. Comparison of the maximum and minimum magnitude 
of gravitational perturbations compared to Saturn's acceleration 



To measure the effect of perturbers on Phoebe's motion, we 
have developed a flexible model. This dynamical model is ad- 
justable, in the sense that it is possible to take into account or not 
the perturbations of the Sun, the eight planets and the Moon, the 
eight major Saturnian satellites and the flatness of Saturn {J 2). 
The masses of the planets are from Folkner et al. ( 2009[> and 
the other parameters of the Saturnian system are from Jacobson 
et al. (2006 ) and are provided in Table [T] 



If a planet is taken into account, it is considered as a point 
mass and its position is from JPL Planetary Ephemeris DE421 
(Folkner et al. 2009 ). In the opposite case, the mass of the planet 
is added to the Sun's mass. Likewise, if a Saturnian satellite is 
taken into account, it is considered as a point mass and its posi- 
tion is given by Lainey et al. ( 2012| l and in the opposite case, its 
mass is added to Saturn's mass. 

For our model, we use a numerical integration of the equa- 
tions of motion using a Gauss-Radau integrator (Everhart 1985 ). 
These equations are written in the Saturnian system barycentre. 
In order to fit the dynamical model, the equations of variation are 
also integrated with the equations of motion (as in |Lainey et al.| 
[2004] l. 

As a large part of astrometric observations is given in ab- 
solute coordinates and not in relative coordinates compared to 
another satellite as for the other main Saturnian satellites, the 
computed position of Phoebe depends on the computed position 
of Saturn. Consequently, the position of Saturn can become a 
source of systematic error in the determination of Phoebe's or- 
bit. 



Name 


Value 




Units 


Phoebe GM 


0.5534 




km 3 .s~ 2 


Mimas GM 


2.5023 




km 3 .s~ 2 


Enceladus GM 


7.2096 




km 3 .s~ 2 


Tethys GM 


41.2097 




km 3 .s~ 2 


Dione GM 


73.1127 




km 3 .s~ 2 


Rhea GM 


153.9416 




km 3 .s~ 2 


Titan GM 


8978.1356 




km .s 


Hyperion GM 


0.3727 




km 3 .s~ 2 


Iapetus GM 


120.5117 




km 3 .s~ 2 


Saturn equatorial radius 


60330 




km 


Saturn/, (xlO 6 ) 


16290.71 






Saturn pole ct P 


40.5955 




degrees 


Saturn pole S P 


83.5381 




degrees 


Saturn polar rate a P 


-0.04229 


dct 


;.century~' 


Saturn polar rate 6p 


-0.00444 


dct 


;.century~' 



In particular, some differences in the position of Saturn be- 
tween several planetary ephemeris exist. In previous models of 
Phoebe, JPL ephemeris DE406 was usually used. Compared to 
the most recent ephemeris, DE421, the difference can reach 0. 15 
arcsec during the 1900-2020 period. In comparison, the differ- 
ence between INPOPlOa ( |Fienga et al.|201 1) and DE421 is less 
than 2 mas during the same period. In these recent ephemerides, 
the orbit of Saturn is determined using Cassini tracking and 
ground-based astrometry. For DE421, the accuracy of Saturn's 
position is tens of km (Folkn eTet al.|2009) . 



3. Observations of Phoebe 

3.1. Ground-based observations 

The ground-based observations of Phoebe come from four dif- 
ferent sources: 

- NSDC: the Natural Satellite Data CenteiQ pr ovides a 
large number of observations of natural satellites ( |Arlot & 
Emelyanov 2009 ). For Phoebe, 3404 observations are avail- 
able from 1904 to 201 1 concerning almost 120 different ref- 



erences; 

MPC: the Minor Planet Center provides observations of nat- 
ural outer irregular satellites of the giant planet^] As of 19 
September 2012, 2040 observations of Phoebe from 1898 
to 2012 are available in the database. The old observations 
before 1981, usually provided in another frame (such as 
B1950.0) have been rotated to the J2000. All the observa- 
tions are given in the J2000 reference frame. Most of the 
observations come from CCD (about 90%) and 5 observa- 
tions made on 13 and 14 June 2010 come from the WISE 



spacecraft (see Sect. 3.2 1; 

NOFS: 205 NOFS observations ( |U.S.N.O| |2012) 
from the Flagstaff Observatory were obtained from 
2000 to 2011. All the data are available on the web 
site of the FASTT Planetary Satellite Observations 
http : //www . nofs . navy . mil/data/plansat . html; 



1 NSDC observations of Phoebe are available on this website : http : 
//www . imcce . fr/hosted_sites/ saimirror/bsapooue . htm 



2 http : //www . minorplanetcenter . net/ iau/ECS/ 
MPCAT-OBS/MPCAT-OBS.html 
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- Pickering: Pickering (1908) published 42 observations of 
Phoebe made in Arequipa, the Harvard observatory station 
from 1898 to 1904. The observations are provided in B1875 
reference frame. They represent the oldest observations of 
Phoebe. The time of observation is given without indication 
of seconds, i.e. to the nearest minute. A specific study of 
these observations has been done in order to improve their 
astrometry (Sect.[4~T). 



Many observations appear in several databases. In particu- 
lar, about 1950 appear in both MPC and NSDC. In that context 
and as a general rule, we favour observations in NSDC because 
they are given in their original format (similar to the publica- 
tion) whereas MPC turns the observations into J2000 reference 
system with their own routine. Consequently, only 83 observa- 
tions from MPC (and only in that database) have been selected. 

In the same way, some observations appear in both the 
NSDC and the NOFS databases. In particular, observations from 
|Stone| ( [20001 |200l) and observations from NOFS compiled in 
NDSC (so0023 and so0029) are replaced by NOFS observations 
( U.S.N.O|2 012 ) that use a more recent reduction process and are 
provided with more digits for time and coordinates. This is also 
the case of MPC observations provided by the NSDC database 
(MPC63362, MPC63363, MP C67130, MPC75 144, MPC75145) 
replaced in our compilation by |U.S.N!O| ( |2012) l. 

Finally, there are also several doubles in the NSDC database. 
Some observations appear in several references but with dif- 
ferent digits for time, or coordinates. In this context, we keep 
the observations with more digits. For example, observations 
in MPC52887 a re replaced by those of NSD C (so0022 file). 
Observations in Mulh olland & Shelus| ( |1980| l appear twice in 
the NDSC database (soOOOl and so0019) and we keep those in 
Strugnell & Taylor] | [1990) (soOOOl). 



3.2. Space-based observations 

During its encounter with Saturn in 1981, Voyager 2 performed 
eight imaging observations of Phoebe. These observations from 
17 June to 15 August 1981 were published in Jacobson (1998). 
They are provided in specific coordinates which are pixel and 
line locations in the Voyager camera frame. The Saturn barycen- 
tric position and the velocity of Voyager 2 during the eight 
observation times are also provided in the paper. The process 
of reduction that takes into account the camera pointing, the 
gnomonic projection, the electromagnetic and optical distortion, 
is fully described in Jacobson ( 1998 1. In this paper, we also use 
pixel and line locations in the camera frame as coordinates and 
the positions of Voyager given in Jacobson ( 1998) are also used. 



The Cassini spacecraft entered Saturn's orbit at the end of 
June 2004. Ever since, its ISS Narrow Angle (NAC) and Wide 
Angle (WAC) cameras have sent tens of thousands of images 
of Saturn, its rings, and its satellites. Right before entering 
Saturn's orbit, Cassini performed a Phoebe flyby on 11 June 
2004, where it reached a minimum distance of about 2068 kilo- 
meters. Cassini's ISS NAC kept observing that satellite from 6 
to 12 June 2004, producing hundreds of high-resolution images 
of Phoebe that can be downloaded from Planetary Data System^] 

A total number of 223 images from Phoebe's flyby were 
used for astrometric reduction. The spacecraft's position and 
camera's pointing vector were computed using SPICE librarjQ 



( |Acton|[T996| l. Since there is an error in the camera's point- 
ing vector, its correction was done using the UCAC2 star cat- 
alogue (Zacha nas et al.|2 004). The satellite's position measure- 
ment was done using the ellipse projection of Phoebe's triaxial 
shape (Thoma s|2010| l. Nonetheless, a difficulty was encountered 
for the centre-of-figure measurement since Phoebe has an irreg- 
ular shape and fitting an ellipse is not the best way to determine 
its position. The astrometric model that was used for pointing 
correction and satellite centre-of-figure measurement is well de- 
scribed in Tajeddine et al. (2013 i. The number of detected stars 
varied from 2 to 43 stars per image, depending on the satellite's 
size on the image and camera exposure length. 

Figure [2] shows an example of Phoebe's astrometric reduc- 
tion from the Cassini spacecraft and the resolution reached. 
With an obtained accuracy of a few kilometers, the Cassini ob- 
servations have proven to be critical in Phoebe orbit modelling. 




Fig. 2. Example of a Cassini astrometric reduction of Phoebe. 
UCAC2 catalogue stars are superimposed on the image. 



Finally, the WISE spacecraft (Wide-Field Infrared Survey 
Explorer) also provided five more space-based observations in 
2009. But unlike to Cassini or Voyager 2, this spacecraft was in 
Earth-orbit and not in the vicinity of Saturn system. These ob- 
servations can be treated as ground-based observations. 

3.3. Catalogue of observations of Phoebe 

A total of 3598 observations (3367 ground-based and 231 space- 
based) from 1898 to 2012 have been compiled in a catalogue. For 

d2007ll used 1606 



http : //pds . nasa . gov/| 



http : / /naif . jpl . nasa . gov/naif / 



comparison, in previous studies Emelyanov 
observations from 1904 to 2007 and |Shen et al.| ( |201 j| t used 2994 
observations from 1904 to 2009. 

For the catalogue of observations, we adopted the format 
presented in Desm ars" et al.| ( |2009b| l and we added some infor- 
mation about the catalogue used for astrometric reduction (us- 
ing the MPC flag), the correction applied on observations, and 
the weight used for the fit. For specific observations, we also 
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add the year of epoch or the angular parameters of the Voyager 
camera. The full catalogue can be requested from the authors. 

Figure [3] represents the distribution of Phoebe observations 
and reveals its strong sparsity. Most of the observations have 
been realized since the mid-1990s. Other observations were from 
the early 20th century and there are a few data between these 
two periods. This strong disparity is a problem for determining 
the accurate orbit of Phoebe. 
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Fig. 3. Logarithmic distribution of the Phoebe observations 



4. Improvement of astrometry 

In this section, we try to improve the astrometry of Phoebe in two 
ways. The first is a new reduction of old data, and the second is 
the correction of bias in stellar catalogues. 



4.1. New reduction of old observations 




In some older publications such as Pickering ( 1908 1 and Perrine 
( |1904| l, authors gave the spherical positions of the natural satel- 
lite and its relative positions to reference stars on the photo- 
graphic plates. In their papers, the position of Phoebe is deduced 
from the positions of reference stars of the Cape Photographic 
Durchmusterung (CPD) catalogue. Unfortunatly, the number of 
reference stars of this catalogue on photographic plates is small 
and their positions are inaccurate. As an illustration, Fig. ^rep- 
resents the statistics of the difference in angular separatiorj^] of 
the stars in CPD and modern catalogues used in current reduc- 
tion. For many stars, the difference in position between the older 
and recent catalogues is more than 5 arcsec. The positions of 
the reference stars in the CPD catalogue represent a source of 
systematic errors on Phoebe's positions. 

Nowadays, although these plates are not available, it is pos- 
sible to reduce the positions of Phoebe from data provided by 
these publications with modern and precise astrometric cata- 
logues such as TYCHO and UCAC4 ( |Zacharias eraT^[20T3 i. 
The authors gave the position of Phoebe relative to the position 



3 4 5 6 7 8 9 10 11 12 13 14 15 

Separation(arcsec) 



5 In this context, angular separation is given by s = 
VAff 2 cos 2 S + A(5 2 where Aa and Ad are the difference between 
right ascension and declination given by older and modern catalogues. 



Fig. 4. Distribution of the separation of pairs of stars in the 
Tycho2-CPD catalogues (top) and UCAC4-CPD catalogues 
(bottom). 



of the reference stars in CPD catalogue. With accurate positions 
of reference stars given in the modem catalogues, it is possible 
to deduce the spherical position of Phoebe on these plates. The 
method of reduction consists in identifying the reference stars 
on the plates and applying corrections of proper motion for each 
star at the date of observation. In the identification process, two 
stars in two different catalogues are considered the same star if 
the distance between their position is smaller than 15 arcsec and 
if the difference of their magnitude is less than one. 

We tried to reduce observations from Pickering ( 19081 with 
this method and by using several stellar catalogues (TYCHO, 
USNOB, and UCAC4). The residuals obtained are at the same 
level of accuracy as the positions provided by Pickering (2-3 
arcsec). In fact, the main difficulty of this method is the poor 
knowledge of the proper motion of reference stars. The observa- 
tions were realized more than 100 years ago and consequently 
inaccurate values of proper motion lead to inaccurate positions 
of stars at the date of observation. For example, Fig.[5]represents 
the position of a reference star given by four different catalogues 
at J2000.0 and B 1900.0 epochs using the star proper motions 
from these catalogues. For the J2000 epoch, the positions of the 
star are quite similar for the different catalogues. For B 1900.0, 
the difference in the positions of the star can reach more than 1 
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arcsec, according to the catalogue used. The HIPPARCOS cata- 
logue is expected to give accurate values of proper motions but 
only a small number of stars from this catalogue is available on 
the plate (less than four). 



-3000 - 
-3500 ■ 
-4000 ■ 



O 
LU 
Q 



-6500 ■ 
-7000 ■ 
-7500 ■ 



■ TYCHO 
USNOB 
UCAC4 

■ HIP 



800 1000 1200 1400 1600 1800 2000 2200 2400 2600 

RA (mas) 

Fig. 5. The positions of one reference star in epoch J2000.0 
(lower) and B1900.0 (upper) calculated with the positions 
and the proper motions in four different catalogues: TYCHO, 
USNOB1.0, UCAC4, and HIPPARCOS. 



The inaccurate values of proper motion of stars are currently 
a difficulty in the new reduction of old data. Finally, for the 
Pickering and Perrine observations, we used the positions pro- 
vided in the papers. A possible solution would be given by the 
Gaia space mission and its astrometric stellar catalogue. Gaia 
will provide accurate positions and accurate proper motions of 
stars (Mig nard et al.|2007] >. In that context, the method described 
in the section could be successfully applied. 



4.2. Stellar catalogue correction 

Stellar catalogues used in astrometric reduction are usually a 
source of systematic error. This error depends on the catalogue 
used and on the zone on the celestial sphere. 

Chesley et al. (2010 1 propose a treatment of star cata- 
log biases in asteroid astrometric observations. Considering the 
2MASS stellar catalogue as the reference, they have detected a 
bias in star positions of five stellar catalogues (Tycho-2, UCAC2, 
USNO B1.0, USNO A2.0, USNO A1.0), deduced from the 
residuals of numbered asteroid observations. They also provided 
a method to remove the biases for observations. 

This method can be applied for the first time to a natural 
satellite such as Phoebe because most of its observations are in 
absolute coordinates, which is not the case for inner or major 
satellites (where inter-satellite positions are also used). Table [2] 
provides statistics on the catalogues used to reduce the astromet- 
ric positions of Phoebe. The stellar catalogue can be identified 
for 2915 observations (86.6%). Among these observations, 2625 
(78.0%) use one of the five catalogues studied inlChesley et al 



(2010) and are pertained to the treatment. Figure [6] represents 



the bias correction applied to Phoebe's observations. As many 
observations of Phoebe have been realized in the same zone on 
the celestial sphere, these observations have the same bias cor- 
rection. Consequently, the number of observations for a specific 



bias correction is also indicated in Fig. [6] The bias removal can 
reach about 0.5 arcsec both in right ascension and in declination. 



Table 2. Statistics on catalogues used for the reduction of the 
Phoebe observations. Code is similar to MPC flag. 
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Fig. 6. Bias in right ascension and declination, and number of 
observations concerned by bias correction (see text). 



5. Fitting process 

5.1. Weighting process 

The fitting process consists in the determination of the epoch 
state vector of Phoebe (components of position and velocity 
at a given epoch) that minimizes the residuals (observed posi- 
tion minus computed position). The classic least-squares method 
(LSM) is used for this determination (for more information, see 
for example Desma rs" et al.|200 9a). In the LSM, a weighting ma- 
trix V b s is required and usually considered as a diagonal matrix 
where the diagonal components are ef - 1/tr? where crj, is the 
estimated variance of the observation i. 
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For our process, we first fit the model to all observations by 
considering their weight as identical and cr — 1 arcsec for each 
observation. Then, the root mean square (rms) of the O-C in 
both coordinates (right ascension and declination) is computed 
for each reference and each observatory. Finally, we fit again the 
model by considering these rms depending on reference and ob- 
servatory as their weight. 

5.2. Correction applied to observations 

In the fitting process, observations can be corrected for different 



Table 3. Mean /j and root mean square (rms) of O-C in arcsec 
for the different models. 



bias. In particular, as described in Sect. 4.2 correction of stellar 
catalogue bias can be corrected. Moreover, for prior 1940 obser- 
vations, a correction in right ascension of 0.75 arcseconds should 
be applied as explained in Jacobson (2000). 



5.3. Change of reference frames 

Most of observations of Phoebe are provided in absolute coor- 
dinates (right ascension and declination) but in a different ref- 
erence system (ICRF, J2000, B1950, true equator and equinox 
of the date, mean equator and equinox of the date, mean equa- 
tor and equinox at 1st January of the year of observation, mean 
equator and equinox at 1st January of a specific year, mean equa- 
tor and equinox of specific epoch). 

Because the dynamical model allows the computation of 
the positions of Phoebe to ICRF, we have to transform obser- 
vations in ICRF in order to compare computed and observed 
positions. The classical transformations have been done with 
the SOFA routine ( |IAU SOFA Board|2010). Fo r transformation 
from B1950, the technique from Murray ( 1989) is applied. 

The case of observations from Pickering (1908) published 
in B1875 is particular. The observations are first corrected 
from elliptical aberration, then transformed to mean equator and 
equinox of the date using the precession value of Newcomb 
( Kinoshita| 1 1 975 [ ). Finally, the coordinates are transformed to 
ICRF by using current values of precession (IAU Resolutions 
2006). 



6. Comparison of dynamical models 

The flexible model developed in Sect. [2] allows the definition of 
four different models of the motion of Phoebe: 
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1.5782 


Model 3 


a 


0.0084 


0.4514 


0.6506 


2.5614 




S 


0.0532 


0.4612 


-0.1232 


1.5793 


Model 4 


a 


0.0079 


0.4675 


0.6897 


2.6214 




6 


0.0458 


0.4649 


-0.1162 


1.5813 



total number of 3367 ground-based observations and all Cassini 
data are accepted, representing 223 observations. 

The residuals for the four models are quite similar in partic- 
ular for Models 1 and 2. As expected, the perturbations included 
in the first model but not in the seconc0have no main influence 
on the motion of Phoebe. 

The residuals for ground-based observations represent about 
4000 km in distance whereas the residuals for Cassini data rep- 
resents only 4.7 km and provide an important constraint on the 
motion of Phoebe (see also Sect. [8]). 

With regard to Table [3] the Model 1 has the best residuals, 
but a simple model with only the perturbations of the Sun and 
Jupiter can be enough to compute the motion of Phoebe. Even if 
the dynamical model is enhanced, there is no clear improvement 
in the residuals. Although the astrometry is better (large number 
of observations, correction of bias, addition of space-based 
observations), it remains unreliable. Consequently, regular 
observations of Phoebe are still necessary. The method of re- 
duction of old observations presented in Sect. |4.1| and the future 
Gaia stellar catalogue also appear to be a good opportunity for 
the improvement of the astrometry. Even if the difference in the 
residuals of the four models are not important, the full Model 1 
provides the best residuals and we deal with this model hereafter. 

Table |4] presents the residuals computed with Model 1 with 
and without correction of the bias in the stellar catalogue (see 
Sect. |4.2[ ). The correction of the bias can help to reduce the resid- 
uals, in particular in declination. 



- Model 1 is the most complete version of the flexible model 
with all the perturbations: gravitational perturbations of the 
Sun, of all the planets and the Moon, then the perturbations 
of the eight main satellites of Saturn, and the flatness param- 
eter (j 2 y, 

- Model 2 includes the perturbations of the Sun, Jupiter, 
Uranus and Neptune, then the three main Saturnian satellites 
(Titan, Iapetus, Rhea), and the fa parameter; 

- Model 3 includes only perturbations of the Sun, Jupiter, and 
Titan; 

- Model 4 includes the perturbations of the Sun and Jupiter. 



Table 4. Post-fit residuals (mean \i and root mean square) after 
correction or not of bias in stellar catalogue. Residuals greater 
than 5 arcsec were rejected. 





Uncorrected 


Corrected 


Number 






rms 




rms 


a 


0.0241 


0.4590 


0.0089 


0.4467 


3319/3367 


6 


0.0857 


0.4856 


0.0564 


0.4617 


3319/3367 



The aim of the four models is to determine if a full model 
is required to compute the position of Phoebe or if a simple 
model is enough. In that context, we compute the residuals of 
the Phoebe observations after fitting to the four models. Table [3] 
provides the mean and the root mean square of the O-C for the 
different models. If O-C were greater than 5 arcsec, they were 
rejected. For each case, 3319 observations are accepted among a 



7. Ephemeris of Phoebe 

In Sect. [6] we have shown that Model 1 including the per- 
turbations of planets and satellites, and the flatness parameter 

6 Gravitational perturbations of Dione, Tethys, Earth, Venus, 
Enceladus, Mars, Mimas, Mercury, the Moon, and Hyperion. 
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provides the best residuals. In this context, Model 1 fitted to 
the whole set of observations is now considered our Phoebe 
ephemeris and is called PH12. The statistics of the residuals for 
each reference is presented in Table [5] Figure [7] represents the 
residuals in right ascension and declination for all ground-based 
observations. 




1880 



1900 



1920 1940 1960 1980 2000 2020 
Date 



n r 



n r 




1880 1900 1920 1940 1960 1980 2000 2020 

Date 

Fig. 7. Residuals in right ascension (R.A.) and declination (Dec.) 
of observations 



We have derived from the numerical integration a binary 
file for ephemeris computation using the Chebychev polyno- 
mials available at ftp : //ftp . imcce . fr/pub/ephem/satel/ 
ph!2/ The PH12 ephemeris is available from 1876 to 2022. 



8. Accuracy of the extrapolated position of Phoebe 

The accuracy of the position of Phoebe can be estimated with 
statistical methods. Desmar s" et al.| ([2009a) presented and de- 
veloped some statistical methods for estimating the accuracy 
of satellite ephemerides. Emelyanovj (2010) also estimated the 
ephemeris precision by three different methods for all outer 
satellites including Phoebe. These methods are Monte Carlo pro- 
cesses and consist in modifying the nominal orbit with specific 
assumptions. The modified orbit is then compared to the nominal 
orbit and the process is repeated many times giving a statistical 
estimation of the accuracy. 

In order to estimate the position uncertainty of Phoebe we 
use two of these methods given in Desmars et al. (2009a). The 
first, Monte Carlo using Covariance Matrix (MCCM) consists 
in adding a random noise to observations, whereas the second, 
bootstrap resampling (BR) consists in randomly sampling the 
observations. For both methods we generate K = 200 clones of 
the nominal orbit. 

A measure of the accuracy in distance cr^ can be computed 
with this relation 



cr d (t) 



\ K 



(1) 



k=\ 



where d*(f) is the distance between the position on the nominal 
orbit and the position on orbit k at time t . 



The position uncertainty obtained by the two methods is rep- 
resented in Fig. [8] Both methods give similar results. The accu- 
racy of the Phoebe ephemeris is quite good during the 2000s 
thanks to a lot of ground-based observations and, in particular, 
thanks to Cassini data. During this period, the accuracy is about 
50km in distance (which represents about 0.01 arcsec in angu- 
lar separation) and less than 10km during the Cassini flyby (June 
2004). Going back to the past, the accuracy deteriorates owing to 
the lack of observations and to inaccurate observations. The ac- 
curacy reaches about 1000km (about 0.2 arcsec) in the early 20th 
century. In the near future, the accuracy will remain less than 
100km (about 0.02 arcsec). As a comparison, Emeiyanov (2010) 
provides an accuracy of 0.05-0.06 arsec during 2010-2020, but 
he used only 1606 observations until 2007 and the Cassini data 
were not fitted. 




1880 



1900 



1920 



1940 1960 
Date 



1980 



2000 



2020 



Fig. 8. Position uncertainty of the Phoebe ephemeris in distance 
calculated with Monte Carlo using Covariance Matrix (MCCM) 
and Bootstrap Resampling (BR) as a function of time in years. 



9. Conclusion 

In this paper we have developed a flexible dynamical model of 
Phoebe that includes the perturbations of the Sun, all the plan- 
ets, the Moon, the main Saturnian satellites, and the flatness pa- 
rameter. We used the most recent theories of planets and satel- 
lites. An important step in the improvement of astrometry has 
been realized. We extended the catalogue of available obser- 
vations of Phoebe from 1898 to 2012, then we dealt with the 
treatment of old observations, and the correction of bias in the 
stellar catalogue. Finally, we have produced a new ephemeris of 
Phoebe. The residuals are 0.45 arcsec and the estimated accu- 
racy of Phoebe's position is less than 100km for the 1990-2020 
period. The comparison of the different models also shows that 
a simple model is enough to determine the position of Phoebe, 
even if the more complete model provides the best residuals. 

The main constraint in the development of an accurate 
ephemeris of Phoebe is the astrometry. Consequently, this is the 
point to improve in the future. One way is to regularly observe 
the satellite. Another way is to improve the process of reduction. 
Recently, Peng et al. (2012) dealt with the geometric distortion 
in astrometric images and they successfully applied their method 
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to the observations of Phoebe. Finally, the main progress for the 
Phoebe ephemeris will be achieved with the Gaia space mission. 
Indeed, Gaia will provide an accurate stellar catalogue (with ac- 
curate positions of stars and accurate proper motions). Thus, the 
method described in Sect. |4.1| could be effective and provide ac- 
curate old data. Moreover, CCD frames and old photographic 
plates could be reduced again with this stellar catalogue. 
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Table 5. Statistics of astrometric residuals. Code bib. is the bibliographic code 
used in the catalogue, n is the number of accepted observations for each refer- 
ence, N is the total number of observations of the reference, /j„ and fi s are the 
mean of the residuals for each coordinate in arcsec, rms„ and rmsj are the root 
mean squares for each coordinate in arcsec and, the time-span is also indicated. 
The symbol * indicates that the observations have been corrected for bias. 



Code bib. Reference 



n/N 



to 



Observations in Strugnell & Taylor (1990) 
Chernykh & Chemykh|(|T97l) 27"2" 0.5854 ^TD816 

Mulholland & Shelus| (l980} 6/8 0.4136 -0.0970 

Observations in 



rms a Time-span 



13 

34 

51 

52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
71 

72 
73 
74 

102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
113 
114 
115 
116 
117 
118 
120 
121 
122 
125 
127 
137 
139 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 



|Pickering|(l904b 
|pickermft||l9U4a} 

|mrine|(l9U4"|> 

|Albrech t & Smitil|( |l909] 
|pemne et aT|(|l9U9|i 
|MNRA S68-1908| fl908 
|MNRA Sb9-19U9| |ffD9 
MJN KAS70-191O iHTTU 



Bec-Borsenberger & Rocher 1 79S2 I 

o/i 



|MN RAS / 1 - 1 9 1 1 1 (p?TT 

Barnard (1913 
Barnard (fPjre 



(1922) 



1944) 



1943) 



van Biesbroeck 
Richmond & Ni cholson 
van Biesbroeck l 
Bobone|(l933J 
van B iesbroeck ( 1957 ^ 
van Biesbroeck ( 193!? 
Koemer & Lloyd ( 196 
van Bi esbroeck el "all 
|Veillet&U ebehogne (TWT 
Debehogne ( 1981 1 
|Smith et all] 1 1981} 
Debehogne H. pnvat c omm. 1981 
iDourneau et al. 



976) 



0/1 
5/5 
11/ 11 
10/ 10 
15/ 15 
23/23 
12/ 12 
7/7 
2/2 
2/2 
0/5 
1/1 
8/8 
7/7 
13/14 
8/8 
2/2 
3/4 
21/21 



{1991} 



i Ledovskaya et a l 
[stone & Harris| p 
|stone||2UUU] i 
Derek Jones coram to NSDC 



-0.4413 
-2.3742 
-1.0285 
-0.1002 
0.6177 
0.7487 
1.3857 
0.4634 
-0.9066 

-0.5528 
1.4202 
-0.7151 
-0.2220 
-0.6959 
1.0106 
-1.2365 
0.5240 

0.3358 
-0.2414 
0.1846 
Observations in NSDC 



18/ 18 
5/5 



|Veiga etal.| ( |2000 i 
MFC comm^2001 -2002a) 
MPC comm. (2001 -2002b) 
MPC comm. (2001 -2002c) 
MPC comm. (2001-2002d) 
MPC comm. (2001-2002e) 
Comm. to NSDC (2002) 
MPC comm. (2002) 
Comm. to NSDC (2003a) 
Comm. to NSDC (2003b) 
Comm. to NSDC (2 003c) 
Fienga et al. 1 2002| 



Comm. to NSDC (2003d) 
Comm. to NSDC (2003-2004) 
Comm. to NSDC (200 4) 
|Peng & Zhang|(|2006| 
|tjiao et al.| ( |2UU6| 
MPC 7065 3 WIS E Observations 
|Qiaoetal.|(2011} 
MPC24160 
MPC38967 
MPC40910 
MPC40911 
MPC41261 
MPC41449 
MPC41830 
MPC43757 
MPC44182 
MPC46761 



9/9 
38/38 
19/ 19 

7/7 
60/60 
22/22 

9/9 

6/6 
18/ 18 

3/3 

3/3 

3/3 
22/22 
38/38 
16/16 
162/ 162 

2/2 
34/34 

5/5 
210/210 
115/115 

5/5 

1173/ 1173 
4/4 
8/8 

38/39 

17/ 17 
3/3 

10/ 10 
3/3 
8/8 

11/11 
6/6 



0.0327 
-0.0901 

0.0158 
-0.0385 
-0.0745 
-0.2601 
-0.3174 
-0.0137 

0.0451 
-0.4311 

0.5031 
-0.0223 

0.0162 

0.0143 
-0.0630 

0.0720 

0.0209 
-0.1025 
-0.0698 

0.0269 
-0.0074 
-0.0593 
-0.0725 
-0.6614 
-0.1696 

0.6800 
-0.1276 
-0.1428 
-0.1008 
-0.2477 

0.3832 

0.4054 

0.0458 



-0.1065 
1.3441 
0.0567 
1.2292 
0.8746 
1.1878 
1.4758 
-0.4752 
-0.5416 

0.0497 
-0.0666 
-0.0210 
0.5337 
1.1637 
0.3198 
0.7450 
-0.1129 

-1.5206 
-0.6121 
-0.2148 

-0.3741 
-0.0090 
0.0829 
0.0701 
0.2488 
0.2854 
0.5602 
-0.0085 
0.1962 
-0.1218 
-0.7376 
0.1123 
-0.4346 
0.1562 
0.0440 
-0.1076 
-0.0003 
0.0402 
-0.0109 
-0.0163 
0.0046 
0.4536 
0.0251 
-0.0325 
-0.0770 
-0.1226 
0.0833 
0.3942 
0.1553 
0.3817 
0.8383 
-0.1473 
0.0652 



1.4508 
0.7573 



0.9992 
2.5230 
1.6793 
0.6667 
0.8672 
1.3078 
1.8302 
1.4150 
1.0271 

0.5528 
1.8178 
1.1774 
0.9805 
0.9347 
1.0171 
2.6553 
0.6780 

0.8243 
0.5502 
0.2326 

0.0795 
0.2336 
0.2770 
0.1007 
0.1828 
0.4536 
0.4106 
0.0318 
0.4876 
0.4765 
0.7787 
0.0453 
0.6929 
1.1432 
0.0857 
0.1855 
0.0423 
0.3632 
0.0772 
0.0592 
0.1191 
0.1572 
0.1486 
1.3239 
0.1981 
1.2349 
0.3540 
0.1443 
0.3738 
0.2626 
1.3833 
0.6700 
0.0739 



1.1405 
0.4710 



0.5282 
2.1368 
0.7945 
1.4391 
1.3982 
1.5321 
1.6171 
0.5778 
0.5457 

0.0497 
0.3589 
0.3068 
0.7550 
1.3354 
0.3413 
0.7780 
0.5557 

1.6113 
0.8142 
0.3068 

0.3846 
0.3486 
0.2154 
0.1803 
0.3986 
0.4563 
0.5755 
0.0306 
0.2463 
0.3222 
0.7420 
0.1131 
1.2151 
1.0772 
0.0564 
0.2209 
0.0036 
0.1568 
0.0693 
0.0575 
0.1508 
0.5128 
0.1099 
0.7528 
0.0933 
1.1008 
0.3465 
0.3948 
0.3201 
0.4060 
1.3944 
0.3502 
0.0988 



1968-1968 
1975-1976 



1904- 
1904- 
1904- 
1905- 
1906- 
1907- 
1908- 
1909- 
1910- 
1912- 
1913- 
1922- 
1940- 
1942- 
1952- 
1955- 
1957- 
1960- 
1969- 
1981- 



1904* 

1904* 

1904* 

1905* 

1908* 

1907* 

1908* 

1910* 

1910* 

1913* 

1913* 

1922* 

1940 

1942 

1952 

1955 

1957 

1960 

1969 

1981 



1981- 1981 

1982- 1982 
1989-1989 



1998- 
1998- 
1999- 
1995- 
1995- 
2000- 
2000- 
2001- 
2001- 
2001- 
2002- 
2002- 
2002- 
1998- 
2003- 
1998- 
2003- 
2003- 
2004- 
2003- 
2003- 
2010- 
2005- 
1994- 
1999- 
1995- 
1998- 
2000- 
2000- 
2000- 
2001- 
2001- 
2002- 



1998* 

1999 

1999 

1997 

1997 

2000* 

2000* 

2001* 

2001* 

2001 

2002* 

2002* 

2003* 

2003* 

2003* 

1999* 

2003* 

2003* 

2004* 

2005* 

2004* 

2010 

2008* 

1994 

1999 

1998 

2000 

2000 

2000 

2000 

2001* 

2001* 

2002* 
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Table 5. continued. 



code bib. Reference 



n/N 



Ha 



rms„ 



mis,} time-span 



210 
211 
212 
213 
215 
216 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 
227 
228 
229 
230 
231 
232 
233 
234 
235 
236 
237 
238 
239 
240 
241 
242 
243 
246 
247 
248 
249 
250 
251 
252 
253 
254 
255 
257 
258 
259 
260 
261 
262 
263 
264 
267 
268 
270 
271 
272 

300 
301 
350 

351 

500 



MPC50594 
MPC51181 
MPC51367 
MPC51497 
MPC53171 
MPC53464 
MPC53629 
MPC53948 
MPC54163 
MPC55508 
MPC55975 
MPC55976 
MPC56148 
MPC56608 
MPC56609 
MPC57946 
MPC57947 
MPC58097 
MPC58523 
MPC58764 
MPC59029 
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MPC69202 
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MPC73667 
MPC74025 
MPC74380 
MPC74811 
MPC75400 
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7/7 0.0202 0.0621 0.0927 0.1233 2004-2004* 

3/3 0.0779 0.0451 0.0801 0.0453 2004-2004 

6/6 0.0582 0.0325 0.0659 0.0418 2004-2004 

4/4 0.0246 0.0406 0.0403 0.0412 2004-2004* 

3/3 -0.1010 -0.0407 0.1257 0.0479 1996-1996* 

19/ 19 -0.0555 -0.0195 0.1914 0.1828 2004-2005 

8/8 -0.0015 -0.0514 0.1990 0.1518 2005-2005 

36/36 0.0253 -0.0210 0.2774 0.2015 2004-2005* 

22/22 0.0228 0.0053 0.2475 0.1784 2004-2005* 

5/5 -0.0147 0.8286 0.2316 0.8815 2005-2005* 

6/6 0.3247 -0.0217 0.7387 0.2284 2005-2005* 

72/72 0.0081 0.0544 0.2003 0.2815 2005-2006* 

34/34 0.0053 0.0140 0.2814 0.2085 2006-2006* 

33/33 0.1470 -0.0966 0.2252 0.1823 2006-2006* 

1/ 1 0.2920 -0.2526 0.2920 0.2526 2006-2006* 

7/7 -0.5558 0.6623 0.8614 1.5035 2006-2006* 

3/3 1.5135 1.2146 1.9990 1.2276 2006-2006* 

10/ 10 0.0283 -0.1981 0.2801 0.4356 2006-2006* 

12/ 12 0.0619 -0.0699 0.5181 0.2336 2006-2006* 

34/34 0.2153 0.0840 0.4006 0.6965 2007-2007* 

9/9 -0.0607 0.1404 0.2501 0.5894 2007-2007* 

23/23 0.1357 0.0140 0.2595 0.2105 2007-2007* 

18/ 18 0.1645 -0.0935 0.2798 0.1918 2007-2007* 

17/ 17 0.2167 -0.1131 0.3165 0.2162 2007-2007* 

14/ 14 0.0379 0.0344 0.3345 0.2506 2007-2007* 

1/ 1 1.2399 -0.2102 1.2399 0.2102 2008-2008* 

3/3 0.3235 -0.1693 0.4012 0.2109 2008-2008* 

30/30 -0.0142 0.0457 0.2321 0.2438 2008-2008* 

20/20 0.0054 0.1463 0.3018 0.4611 2008-2008* 

23/23 0.0406 -0.0538 0.2190 0.2306 2008-2008* 

8/8 0.2757 0.1044 0.4027 0.1738 2008-2008 

2/2 0.4133 -0.0853 0.4201 0.0967 2008-2008 

14/ 14 0.3811 -0.2029 0.4459 0.2786 2008-2008 

3/3 0.2071 0.2070 0.3299 0.2858 2008-2008 

4/4 -0.0818 0.3491 0.4814 0.5066 2008-2008* 

6/6 0.0161 0.5311 0.2223 0.6249 2008-2008* 

9/9 0.0818 0.0208 0.1599 0.0503 2009-2009* 

4/4 -0.5536 0.2077 0.6222 0.2978 2009-2009 

15/ 15 0.1077 -0.0344 0.2934 0.2110 2009-2009 

16/ 16 -0.1531 0.3044 0.5297 0.8047 2009-2009* 

2/2 -0.5985 0.0678 0.6111 0.3130 2009-2009 

6/6 0.0406 0.0375 0.2085 0.1453 2009-2009* 

4/4 0.1597 -0.1154 0.2672 0.2136 2009-2009 

2/2 -0.0147 -0.1053 0.0530 0.1070 2010-2010* 

6/6 -0.0511 0.0803 0.2327 0.1850 2010-2010* 

30/30 0.0601 0.3838 0.3508 0.8498 2010-2010* 

21/21 0.6952 0.9287 0.9956 1.3787 2010-2010* 

4/4 0.0733 0.1649 0.1064 0.1968 2010-2010* 

1/ 1 0.0851 0.1336 0.0851 0.1336 2010-2010* 

2/2 -0.0898 0.0278 0.1128 0.0305 2011-2011* 

21/21 0.0002 -0.0429 0.1241 0.1678 2011-2011* 

39/39 0.2880 0.1457 0.5760 0.5211 2011-2011* 

29/29 0.1555 0.0254 0.3181 0.1771 2011-2011* 

3/3 0.2502 0.0030 0.2608 0.0727 2011-2011* 

9/9 0.3081 -0.0401 0.3166 0.1026 2011-2011* 

6/6 0.0941 -0.2012 0.1434 0.2334 2011-2011* 
Other observations 

6/42 0.2576 0.9413 1.9223 3.0575 1898-1902* 

205/205 0.0201 0.0420 0.2786 0.3851 2000-2011* 

8/8 -6.6018 -1.1079 6.6194 1.2836 1981-1981 

223/223 0.6304 -0.1291 2.5353 1.5782 2004-2004 

83/83 0.1518 0.1092 0.4838 0.3341 1907-2012* 



